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1. Introduction 



The interest in the heavy quark spectroscopy has been revived due to the 
discovery of new heavy hadrons, containing the charmed quarks, over the 
last few years [1, 2]. In 2003 a new resonance named A (3872) was reported 
by the Belle Collaboration in the invariant mass distribution of J/ipn+TT^ 
mesons produced in B ± — > A ± A(3872) — > K ± J/ip7r + 7r~ decays. It appeared 
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as a narrow peak with a mass 3871.2 ± 0.5 MeV and a width T < 2.3 MeV 
[3]. This state was confirmed by BaBar [4], CDF [5] and DO Collaborations 
[6]. 

In July 2005 Belle claimed the observation of a new charmonium reso- 
nance named X(3940) with a mass of 3943 ± 6 ± 6 MeV and a total width 
of less than 52 MeV [7]. 

In the early 80's, Gelmini [8] studied the S'-wave ccnn states using the one- 
gluon-exchange potential and the virtual annihilation of color pairs, obtaining 
some candidates that could lie below any of the dissociation channels. Chao 
[9] explored the decay, hadronic production, production e + e~ annihilation, 
and photoproduction of various types of ccnn states, using the quark-gluon 
model proposed by Chan and Hogaasen [10]. After the discovery of the 
X(3872) the question about the possible existence of ccnn bound states was 
posed again. Maiani et al. [11] constructed a model of the X(3872) in terms 
of diquark-antidiquark degrees of freedom. Using the X(3872) as input they 
predict other ccnn states with quantum numbers ++ , l + ~ and 2 ++ . Ebert 
et al. [12] addressed heavy tetraquarks with hidden charm in a diquark- 
antidiquark relativistic quark model, concluding that the X(3872) could be 
identified with the 1 ++ neutral charmed tetraquark. 

The consideration of relativistic effects in the composite systems is suffi- 
ciently important when the quark structure of the hadrons is studied [13-21]. 
The dynamical variables (form factors, scattering amplitudes) of compos- 
ite particles can be expressed in terms of the Bethe-Salpeter equations or 
quasipotentials. The form factors of the composite particles were considered 
by a number of authors, who have in particular applied a ladder approxi- 
mation for the Bethe-Salpeter equation [22], ideas of conformal invariance 
[23] , a number of results was obtained in the framework of three-dimensional 
formalisms [24] . It seems that an application of the dispersion integrals over 
the masses of the composite particles may be sufficiently convenient to the 
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description of the relativistic effects in the composite systems. On the one 
hand, the dispersion relation technique is relativistically invariant one and 
it is not determined with a consideration of any specific coordinate system. 
On the other hand, there is no problem of additional states arising, because 
contributions of intermediate states are controlled in the dispersion relations. 
The dispersion relation technique allows to determine the form factors of the 
composite particles [25]. 

In papers [26, 27] the nucleon form factors were calculated in the dis- 
persion relation technique. In paper [26] the relativistic generalization of 
the Faddeev equations was constructed in the form of dispersion relations in 
the pair energy of two interacting particles and the integral equations were 

obtained for the three-particle amplitudes of S'-wave baryons: for the octet 

1+3+ 
J p = — and the decuplet J p = — . The approximate solution of the rela- 
tivistic three-particle problem using the method based on the extraction of 
the leading singularities of the scattering amplitudes about sn~ = 4m 2 was 
proposed. The three-quark amplitudes given in Refs. [26, 28] could be used 
for the calculation of electromagnetic nucleon form factors at small and inter- 
mediate momentum transfers [27]. Using these results the proposed approach 
was generalized to the case of five particles in the recent paper [29]. 

In the present paper the computational scheme of the electromagnetic 
form factors of the tetraquarks, consisted of four particles, in the infinite 
momentum frame is given. 

Section II is devoted to the calculation of electromagnetic form factors 
of N particle system and a special case N = 4 in the infinite momentum 
frame. The calculation results of electric form factors of the tetraquark states 
X(3872) and X(3940) are given in Section III. The last section is devoted to 
our discussion and conclusion. 
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2. The calculation of electromagnetic form factors of tetraquarks 
in the infinite momentum frame 



The approach used in [29] for considering of form factors of exotic baryons 
with isospin 1 = 5/2 and based on the transition from Feynman amplitude 
to the dispersion integration over the masses of the composite particles, may 
be extended for multiquark systems containing N quarks. 

Let us consider the electric form factor of a system of N particles, shown 
in Fig. la. The momentum of the system is treated to be large: P z — > oo, the 
momenta P = ki + k 2 + ... + k^ and P' = P + q correspond to the initial 
and final momenta of the system. Let us assume P = (Pq, P± = 0, P z ) and 
P' = (Pq, P'j_, P z ), where P 2 = s, P' 2 = s'. Then we have some conserva- 
tion laws for the input momenta: 



N 



i=l 

N 

P z -J2k iz = P z (l-J2 Xi )=0, 
1=1 1=1 

N N 1 / N 2 \ 

1=1 i=l z \ 1=1 1 J 



=1 

N N 



i=l i=l 

«■! = rn 2 + k 2 ± , Xl = ^, i = 1,2, ...,7V. (1) 
By analogy for the output momenta : 

N 



k' 1± + ^k J± -q ± = 0, 

N N 

Pz ~ k\ z — ki z = P z [z — x x — = 0, 

i=2 i=2 
i=2 i=2 2 \ 1 i=2 1 J 

A = % m' 2 ± = mj + K\. (2) 
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P' g> _|_ g _ g 2 

It is introduced in (1) and (2) qj_ = P'j_ and z = = . The 

form factor of the iV-quark system can be obtained with the help of the 
double dispersion integral [25]: 

(3) 



F(q 2 ) = J 



ds ds' disc s disc s 'F(s, s', q 2 ) 
(s - M 2 ){s' - M 2 ) ' 



N 

(E m i) 2 

i=l 



disc s disc s 'F(s, s' , q 2 ) = GG' J dp(P, P', hi, k 2 , A;jv-i) (4) 

The invariant phase space dp(P, P', k±, k 2 , fcjv-i) , which enters in the dou- 
ble dispersion integral, has the form [25]: 

dp(P, P', k ± , k 2 , fcjv-i) = d& N \P, h, k 2 , k N )xd& N \P', k[, k' 2 , k' N ) X 



N 



x [J(27r) 3 2fc i0 5 3 (ki - kj), 



(5) 



1=2 

where the iV-particle phase space is introduced: 

N N 

d*W(P, h, k 2 , k N ) = (2n)^(P -E**)II ( 27 r)32(^o) 2 ' 



d 3 k, 



After the transformation we have: 



dp(P, P',h, k 2 , ...,k N -i) = 



( 



2 Ar - 1 (27r; 



3N-5 



N ~ 1 A 1 

n -*u x — \ 



1=1 



Xi 
\ 



z — 1 + X\ 



■X 



For the diquark-spectator (Fig. lb) the invariant phase space takes the fol- 
lowing form similar to (6): 
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(7) 



where the phase space of the diquark is determined by In-i,n- In such a way 
the expressions for any number n (n < [^p-]) of diquark-spectators may be 
obtained. 

Let us consider the latest 5-function in (6). Assuming that: 
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it is obvious that: 
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i± 
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N± 



--2 Xi 
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= 5 



s' + ql 



1 - E x i 1 
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m 



i_i_ 



s + 



Z Z — 1 + X\ Xi 

So we get an expression, which doesn't depend on N. Hence the integration 
over s' is carried out exactly in the same way as in the previous paper [29]. 
The expressions for s and s, which are obviously depend on the number of 
particles N, have the following form: 
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N-2 j-lJV-2 , 

n-2 2 rn%_ 1N + E fcu + 2 E E a/ fci±*y± cos (0; 

2-^1 rr ■ N-2 

i- $>* 

i=i 

JV-1 j-lJV-l 



+ ^ (g) 
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Having integrated the (^-functions, we obtain the following expression for N- 
particle contribution to the form factor (3): 

2 K y N-i i jv-i y iv-i 

f « = 2 2 (j v-i) (27r) 3(7v-i) y n ^ J U dx j u d & x 

4-1 1-1 4-1 

1 bX+l 0(A, - s)0(A, - g) 

1=1 

the following one is for the case with one diquark-spectator: 



k± N-2 \ N-2 2l l N-2 



X 



Fi = 2 2 {N ^^(n-d y U dk i± y n^y 

4-1 4-1 l ~ l 

1 &A + 1 fl(A. - s)6(A s - sQ 
^2 b + \ f (s — M 2 )(V — M 2 ) 1 J 

i=l 

6 = * 1 + ^, / = 6 2 - 4 ^ C ° s2(0l) , 



-^^ 2 -/)(l-(J)/) 



A = * j * , s> = s + q 2 (l + 2A); (11) 
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Akf ± cos 2 (0i) 



A = ! ^ ^— , 5' = 5 + g 2 (l + 2A). (12) 
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To find the tetraquark form factor one needs to put N = A and to account 
the interaction of each light quark with the external electromagnetic field us- 
ing the form factor of nonstrange quarks f q (q 2 ) [30]. For the electromagnetic 
form factor of tetraquark in the case of the normalization G E (0) = 1 we 
obtain: 



U {q } F*(0) f q (0) J 6 (0) + J 9 (0) ' 1 6} 



where: A k . i 

2 * 2 „ 2 7T 2 



X 

X2) 



b\ + l 2 0(A 8 -s)6(A 8 - s') 
X b + \ f 3 {s-M 2 ){s> -M 2 )' 

A k , „ 1 



■W) /II /II «fe f" II * ,(!-„)„(!'„)„(!- ' X 



£3) 



6A + 1 2 g(A.-5)g(A a -gQ 

b + Xf 1 (s — M 2 )(s' — M 2 ) ' 1 J 

A n (n = 1,3) determine relative contributions of the subamplitudes in the 
total amplitude of the tetraquark [31]. 



3. Calculation results 



The electric form factor of tetraquark is the sum of two terms (13). The 
phase space of the spectator contributes to the first term 7 34 = 9.11 GeV 2 . 
The vertex functions G and G' are taken in the middle point of the physical 
region. The mass of the quark u is equal to m u = 0.385 GeV, and the 
mass of the c quark is m c = 1.586 GeV. The cutoff parameter over the pair 
energy A = 10 is obtained in [31]. It is possible to calculate the dimensional 
cutoff parameters over the total energy and the transvers momentum A s = 75 
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GeV 2 , A k 2^ = 1.5m 2 (m 8 = m u , m c ) respectively. It is necessary to account 
that 7-quantum interacts only with the light dressed quark, which has own 
form factor [30]: for u quark f q (q 2 ) =exp(a q q 2 ), a q = 0.33 GeV -2 . We 
can use (13) for the numerical calculation of the tetraquark form factor. 
It should be noted, that the calculation has not any new parameters as 
compared to the calculation of the meson and tetraquark mass spectrum. 
In the previous paper [29] proton charge radius was calculated. It turned 
out to be R p = 0.44 fm, that is almost a factor of two smaller than the 
experimental value R p e xp = 0.706 fm [32]. It is usually for the quark models 
with the one-gluon input interaction [33, 34], when only the presence of the 
new parameters or the introduction of an additional interaction allows to 
achieve a good agreement with the experiment [35, 36]. 

The behaviour of the electromagnetic form factor of the tetraquark X (3872) 
with the mass M = 3872 MeV is shown in Fig. 2. The calculations were car- 
ried out for two tetraquarks with the masses 3872 MeV and 3940 MeV. The 
results turned out to be equal: the charge radius of the X(3872) and X(3940) 
Rtetra = 0.50 fm.The charge radius was found to be approximately equal to 
the charge radius of the proton (R p = 0.44 fm) [29]. It can be concluded that 
tetraquarks with charm are more compact systems than ordinary baryons. 

4. Conclusion 

The method applied in the present work based on the transition from the 
Feynman amplitude to the dispersion integration over the masses of the com- 
posite particles was extended to the system of N quarks for the multiquark 
states and was applied in special case of four particle systems - tetraquarks. 
The absence of any new parameters introduced in the model for the compu- 
tation of the tetraquark form factors is an advantage of this method. 
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Figure captures 

Fig.l Triangle diagrams, which determine the form factors of tetraquarks. 

Fig. 2 The electromagnetic form factor of the tetraquark X(3872) with 
mass M=3872 MeV. 
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